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Plasma Instabilities in an Anisotropically Expanding Geometry 
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We study (3+l)D kinetic (Boltzmann-Vlasov) equations for relativistic plasma particles in a one- 
dimensionally expanding geometry (a special Kasner-type universe) motivated by ultrarelativistic 
heavy-ion collisions. We set up local equations in terms of Yang-Mills potentials and auxiliary 
fields that allow simulations of hard-(expanding)-loop (HEL) dynamics on a lattice. We determine 
numerically the evolution of plasma instabilities in the linear (Abelian) regime and also derive their 
late-time behavior analytically, which is consistent with recent numerical results on the evolution of 
the so-called melting color-glass condensate. We also find a significant delay in the onset of growth 
of plasma instabilities which are triggered by small rapidity fluctuations, even when the initial state 
is highly anisotropic. 



Plasma instabilities have recently been suggested to 
play a major role in the equilibration of matter created 
by an ultrarelativistic heavy-ion collision, e.g. at the Rel- 
ativistic Heavy Ion Collider (RHIC) or the Large Hadron 



Collider (LHC) 0, S, 1 11 [ 
lision, saturation scenarios M, 



Shortly after such a col- 
indicate that the typical 



momentum of a particle in the local plasma rest frame 
is much larger transverse to the collision axis than par- 
allel to it. This momentum anisotropy inevitably leads 
to a so-called Weibel Q (or filamentation) plasma in- 
stability that manifests itself by rapidly growing trans- 
verse magnetic fields. If large enough, a transverse mag- 
netic field bends the trajectories of particles out of the 
transverse plane, thus making the system more isotropic. 
Plasma instabilities are therefore a prime candidate for 
causing rapid isotropization of a quark-gluon plasma, es- 
pecially since they act on a time-scale that is paramet- 
rically shorter than that of scatterings by at least one 
power of the strong coupling constant g. 

There are, however, some caveats: for instance, numer- 
ical simulations 0, @] of non-Abelian Yang-Mills equa- 
tions in the (stationary) hard-loop approximation [lOj 
have shown that the initially exponential growth of gauge 
fields slows down to a weak linear growth when self- 
interactions of the unstable modes become important and 
energy in unstable modes cascades to stable modes of 
higher momentum Interestingly, this may lead to 

the generation of anomalously low viscosity [12j. 

Maybe more importantly, the matter created in a 
heavy-ion collision is believed to escape relatively unim- 
peded in the longitudinal direction (the direction of the 
collision axis). This effectively one-dimensional expan- 
sion decreases the density of hard particles 13, 3] and 
thus attenuates the growth rate of plasma instabilities. 
On the other hand, the expansion increases the de- 
gree of anisotropy, thus making more and more higher- 
momentum modes unstable. While numerical simula- 
tions of classical Yang-Mills dynamics have provided 
some qualitative understanding of the counterplay of 
plasma instabilities and expansion [l5| , an analysis based 



on the hard-loop approximation is desirable, in partic- 
ular to address systematically the fate of non-Abelian 
plasma instabilities and of the associated energy cascade. 
The purpose of this Letter is to develop the basis for 
such a treatment. We generalize the anisotropic hard- 
loop effective theory of Refs. 0, Q to the case of a dy- 
namical, boost-invariantly expanding background, thus 
preparing the ground for corresponding lattice simula- 
tions. We work out explicitly the already rather nontriv- 
ial dynamics in the linear (Abelian) regime (which thus 
are in principle of interest also to ultrarelativistic conven- 
tional plasma physics) and discuss possible implications 
for ultrarelativistic heavy-ion collisions. 

Ignoring the effects of collisions, the dynamics of col- 
lective modes in a non-Abelian plasma is determined by 
the gauge covariant Boltzmann-Vlasov equations 



P 



■ D 6 f a (p, x, t) = gPuFrdl p) f (p, x, t), 



coupled to the Yang-Mills equations 
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-5f a (p,x, t). 
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Here Fjf" = d^A v a - d v 'A£ + gf abc A£A» is the field 
strength tensor in the adjoint representation, D£ c = 
Sued 11 + gfabcAfr is the gauge-covariant derivative, A 11 is 
the gauge-potential. /o(p,x, t) is the (suitably normal- 
ized) color-neutral background distribution of hard par- 
ticles and 6f a (jp,x, t) its colored fluctuations. Eq.|T|) re- 
quires that the background /o satisfies p- d /o(p, x, t) = 
(which is trivially fulfilled when /q is space-time indepen- 
dent). Assuming that the matter created after a heavy- 
ion collision expands longitudinally in a boost-invariant 
way and is sufficiently large in the transverse direction 
such that transverse gradients are small, we take 

fo(p,x) = fo(p±,p z ,z,t) = fo(p±,tp z - zp°) (3) 

which also satisfies p ■ d /o(p, x, t) = 0. 

To describe the dynamics of fluctuations around a 
boost-invariant background, convenient coordinates are 
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proper time r = \/t 2 — z 2 and space-time rapidity 77 = 
arctanh|. The metric in these coordinates is g T r= L 
gij= —Sij, g V ri= —t 2 1 i=l,2, thus corresponding to a 
one-dimensionally expanding space-time geometry. 

Transforming the set of equations (|1I2|) to the coordi- 
nates x a — (r, , 77) (denoted collectively by Greek let- 
ters from the beginning of the alphabet) we find 



P-D6f a \ p , 



gp a F^ p) f (p ± ,p v ), 



-D a {tF' 

T 



a/3\ _ 



J 2 



d 2 pj_dy 
(2tt) 3 



(4) 
(5) 



where the derivative of 5f a has to be taken at fixed p M 
as opposed to fixed p a . Here we have introduced also 
momentum rapidity y = atanh(p°/p z ), such that p T — 
\p ± \ cosh (y — 77), p n — Ip^lr^ 1 sinh (y — rj). 

We assume that at the hypersurface r = T- lao the back- 
ground /o is isotropic and choose the specific model 



fo{p,x) = f iso (^Jp\ + P 2 /t? so ^ 



(6) 



which corresponds to increasingly oblate momentum 
space anisotropy at r > t iso (but prolate anisotropy for 
t < Ti SO ). However, in what follows we shall start the 
time evolution at nonzero proper time To, allowing for 
the fact that a plasma description will not make sense 
at arbitrarily small times, and we shall mostly consider 
the situation that the initial momentum distribution is 
highly oblate at To, i.e. To 3> t; so . The distribution func- 
tion fo has the same form as the one used in Refs. @, Hj], 
but the anisotropy parameter £ therein is now space-time 
dependent according to £(r) = (T/T iso ) 2 — f . 

Since p ■ d [d^L\fo(px,p v )] = 0, we can solve Eq. (@]) by 
introducing auxiliary fields W£(r, x l , 77; v l , y) in 

Sr(x;p) = -gW^T,x\ m v\y)d^ ) fo(px,Pr 1 ) (7) 
that obey 



v ■ DW a (T,x\ri;v\y) = v p F a0 



(8) 



where v a = t^-j = (coshfa — 77), cos cf>, sin <ft, smh fa 21 ). 
At any given space-time point the fields W a depend only 
on the velocity of the hard particles and not on their 
momentum scale, and thus directly generalize the auxil- 
iary fields W fl (x; v) of the hard- loop formalism in a static 
background distribution fo [13, flif . 

For fo of Eq. (J6j) , the induced current j 13 takes the form 
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' i, 



2tv 
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dy v a 



v J Wi 



W r . 



(9) 

ff 2 /o°° ^5-/iso(p) is the Deb y e mass of 

the isotropic case. A short calculation confirms that the 
current is covariantly conserved, D a rj a = 0. 



where m 2 D 



The Eqs. ©, © together with the Yang-Mills equa- 
tions (O can be simulated on a 3-dimensional lattice by 
discretizing space-time and introducing lattice links in 
a standard way [13] and by also discretizing the resid- 
ual momentum variables <f> and y in order to have a fi- 
nite number of W fields. This directly generalizes the 
discretized hard-loop effective equations of motions of 
Ref. [i[ to what may be called the hard-expanding-loop 
(HEL) case. As in the stationary case, momentum dis- 
cretizations that respect reflection invariance (now with 
respect to <p and y — rf) automatically ensure covariant 
current conservation, D a Tj a — 0. 

In this Letter, we initiate this program by studying the 
onset of plasma instabilities and their evolution in the lin- 
ear regime, where their non-Abelian self-interactions are 
still negligible and where we can avoid a discretization of 
cj> and y by solving the equations of motions of the auxil- 
iary fields W a (T, x\ 77; (j>, y). In a stationary plasma with 
oblate momentum-space anisotropy, the most unstable 
modes have wave vectors along the longitudinal direc- 
tion. A particularly interesting case can thus be studied 
by neglecting the transverse dynamics (diA a — 0). Lin- 
earizing in the gauge potentials we have 

[T- l d T Td T - T- 2 d 2 ] A\ T , 77) = f, TdrT^drA,, = j„ 

(10) 

in the gauge A T = 0. Eq. ([8]) can be solved by the method 
of characteristics which gives 



W a (T,7r,(/),y) 



dr 



V F a B I , , ,. 
/ ! It' ,r}{r') 



cosh(?/ — ^(t')) 
y — t](t') — asinh^— sinh (y — r/)j , (II) 

where t](t') is the solution along the characteristic. 
Within our approximation we can thus proceed to eval- 
uate the dtp integral in Eq.®, finding 

/ „ \ -2 
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dy l 



T 2 
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dr' 



t' 'J t^ cosh 77' 



D 



dy v n v ri 



dr' d r i A^fr' , 77'), (12) 



where 77' = 17 (r') and 77' = ?7(t') — y. 

Introducing a Fourier transform in space-time rapidity, 

dv 



A\t,t}) = 



2tt 



exp(ivr])A l (t, v) 



(13) 



and choosing ^(tq) = for simplicity, we find 
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m 2 D f dy sinh y 



dr'A^r', v)d T ,e lu \- y ~™ in ^ sinhy )} X. (14) 
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Below we solve the integro-differential Eqs. (|10p . 
(H3J) numerically. The late time behavior t » To, 
however, may be studied analytically by expanding 
exp \iv (y — asinh(p- sinhy))] in Eq. (|14p around r' = r 
and subsequently acting with d 2 r 2 on Eqs. fTU)) . In this 
limit, the integro-differential equations turn into ordinary 
differential equations for each mode v, 

[d 2 rd T Td T + v 2 d 2 r + fid 2 T - h^t- 1 ] A 1 (t, v) ~ 0, (15) 

[d T T~ l d T + 2/iT" 2 ] 2, (r, v) ~ 0, (16) 

where fi = ^m 2 D TTTi so . From Eq. (|16p we find for the late- 
time behavior of longitudinal fields A z = t~ 1 A v coshry 



_1 A(r, v) 



(17) 




1 5 10 20 30 50 



3G0 



FIG. 1: Numerical solution for individual magnetic modes 
B 1 (t,v) normalized to E 1 (t ,v) = r _1 II ! (ro, v) (thick dark 
lines) for To/n so = 100 and c = 0.5 in a log-square root plot. 
Thin light lines give the analytical result of Eq. (|21[l with nor- 
malization to match the numerical result at late times. Note 



that r, 



o 



Q s ~ 1 and 3 GeV for RHIC and LHC, resp. 



where J n (x) and Y n (x) are Bessel functions of the first 
and second kind, respectively, and Ci,2 are constants. The 
asymptotic behavior of these functions is oscillatory with 
amplitude ~ x~ x / 2 , so t _1 A n only has stable modes. 

For very infrared modes i/<1, the terms proportional 
to v in Ea.(|15p can be neglected and we find for A 1 

A\t, ~ ci Jo (2^FF) + c 2 Y (2V/F) , (18) 

which again correspond to stable oscillatory solutions, 
but whose frequency is a factor \J~2 smaller than that 
of A v . This is consistent with the analytic results of 
Refs. [2, 4] on nonexpanding anisotropic plasmas, where 
the ratio of the longitudinal and the transverse plasma 
frequency was also found to approach v2 in the limit of 
infinite anisotropy parameter £. 

For high- momentum modes v 3> 1, however, only the 
terms proportional to v in Eq. lllSp matter, and we find 

2 l (r, ~ Cl yf?h (2^F) + czyfiFKi , (19) 

where /, K are modified Bessel functions with asymptotic 
behavior I n (x) ~ cx p( x ) K n (x) ~ cx j^ x ^ . Clearly, the 

first of these solutions corresponds to a rapidly growing 
mode which leads us to expect that the large v modes of 
A % will be the dominant modes at sufficiently late times 
with a behavior of 



A\t) ~ t 1 / 4 exp (2^/17) . 



(20) 



This behavior has qualitatively been found by numeri 
cal simulations of the melting color-glass condensate 
For moderate momenta v ~ 1 the solutions to Eq.JT 



are more complicated and can be given in terms of gen- 
eralized hypergeometric functions 2-P3 and a Meijer G- 
function. The dominant contribution turns out to be 

2 l (r, v) ~ t 2 F 3 (^, 2±£; 2, 2 - iv, 2 + iv; -/it) (21) 



with s = vl + 4^ 2 , which interpolates between the sim- 
ple cases 1/ < 1 (Eq. (fig])) and 1/ > 1, Eq. ([201 . 

Information on the behavior of A 1 at very early times 
can be gained by studying Eq. (fl"4"]) for v ^> t/t; so . Ex- 
panding around the stationary point y = 0, we find the 
oscillatory behavior A % (t, v) ~ c\T %v . 

In order to investigate the onset of plasma insta- 
bilities we have solved the non-approximated integro- 
differential Eqs. (|10p . (|14p numerically. To this end we in- 
troduce n*(r, v) = Td T A z (r, v) which obeys d T IP(T, v) = 
—v 2 t~ 1 A 1 (t,v) + rj l , and apply a leap-frog algorithm 
to solve the coupled equations after discretizing both the 
variable r and t' of the memory integral in f , with initial 
conditions A 1 (to, v) — 0, II 1 (to, v) ^ for a given v. 

In order to fix our dimensionful parameters in a way 
that makes contact with heavy-ion physics, we adopt 
the saturation scenario 0, Q and assume that at To ~ 
Qj 1 we have an initial hard-gluon density [14J ti(to) = 
2cQl/(3ir 2 a s Q s T ), where Q s is the saturation scale. 
We consider two cases for the gluon liberation factor c: 



c ~ 0.5 according to numerical simulations of Ref. [17 1 
and c = 2 In 2 according to an approximate analyti- 
cal calculation of Ref. [18]. Assuming further a de- 
formed thermal distribution /o with Ti SO < to and trans- 
verse temperature T = Q s /d with d^ 1 = 0.47 taken 
from Ref. Q, we obtain for a purely gluonic system 
tx/Q a = cd7r 2 /[48C(3)] « 0.182 for c = 0.5 and « 0.505 
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TABLE I: Approximate values of the proper time rio where 
the first of the modes A 1 (r, v) has grown by a factor 10. 



To/r lS o 




+ 1 


1 


10 


100 


1000 


TlO/Vo 


(c 


= 0.5) 


260 


60 


50 


49 


tio/to 


(c 


= 2 In 2) 


95 


25 


21 


20 



for c = 2 In 2. The amount of anisotropy at time tq is 
determined by the ratio tq/t iso . 

In Fig. [1] we display our numerical results for the mag- 
netic field strength B % (t, v) = vt^ 1 A 1 (t, v) with c = 0.5 
for v = 3, 10, 30 and high initial anisotropy, tq/t iso = 100 
(i.e., £(to) ~ 10 4 ). For late times we observe a nearly per- 
fect agreement with the analytical estimate of Eq . (f2"Tj) . 
which however does not contain information on the am- 
plitudes of the unstable modes compared to their ini- 
tial values. In the numerical evaluation we find that the 
modes with larger longitudinal momentum v 1 which have 
larger growth rates at late times, typically start to grow 
also only at larger proper times. This causes a certain 
delay for the onset of the instability which can be mea- 
sured, e.g., by the time t\q when the first of the modes 
Ai has grown in amplitude by a factor of 10 relative to 
its first maximum. For r > t\q the dynamics of the 
collective modes is clearly dominated by the unstable 
modes. For times r up to about half of rio we observe 
a decrease of the energy density carried by the unstable 
modes due to the expansion, but increase despite contin- 
ued expansion thereafter. In Table H] we list the values of 
Tio that we found by scanning through v with various ini- 
tial anisotropies. Ideal-hydrodynamic fits to experimen- 
tal data at RHIC indicate a life-time of a quark-gluon 
plasma of less than 5 fm/c ~ 30 To [l9(. The values ob- 
tained in Table U turn out to be too large to suggest an 
important role of plasma instabilities in RHIC experi- 
ments if the gluon liberation factor c ~ 0.5 as obtained 
in Ref. [3, 17|. However, these unstable modes could 
contribute to fast isotropization of a quark-gluon plasma 
if the value of c is much larger than those currently con- 
sidered, or if the viscosity of the plasma is significant (in 
which case the life-time is increased 20]). Finally, our 
results suggest that even for c ~ 0.5, plasma instabilities 
will be an important phenomenon at the LHC, where 
plasma life-times might exceed 7 fm/c~ 100 To 21 1. 

To conclude, in this Letter we have laid the basis for 
numerical simulations of non- Abelian dynamics in heavy- 
ion collisions by generalizing the anisotropic but station- 
ary hard-loop effective theory to the longitudinally ex- 
panding case. This opens the possibility of studying 
the fate of non- Abelian plasma instabilities in the highly 
nonlinear regime with expansion included. Moreover, we 
have determined the evolution of Weibel instabilities in 
a longitudinally expanding plasma in the regime where 
non-Abelian self-interactions of the unstable modes are 



negligible, finding that they start out oscillatory and later 
grow fast with an asymptotic behavior that we could de- 
termine analytically. Numerically, we have also been able 
to quantify the onset of this instability. Matching our 
dimensionful parameters to those of the saturation sce- 
nario we find that plasma instabilities overcome the ef- 
fects of expansion at or after r ~ 20Q7 1 (which is roughly 
consistent with Ref. [l5| ) . Unless initial parton densi- 
ties are significantly higher than assumed here, only the 
prospected LHC experiments seem to offer large enough 
Q s and plasma life-times to generate strong quark-gluon- 
plasma instabilities from small rapidity fluctuations. The 
dynamics of strong initial fluctuations can only be deter- 
mined by full nonlinear studies; this is work in progress. 

We are grateful to Mike Strickland for collaboration 
and stimulating discussions. PR was supported by 
BMBF 06BI102. 
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